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We apply a Hartree-Fock approximation to a two-orbital model proposed for Fe pnictide 
superconductors. It is found that the antiferromagnetic (AFM) order with the ordering vector 
Q — (tt, 0) is realized. The AFM order appears simultaneously with ferro-orbital order, the 
latter leads to a secondary lattice distortion. We also investigate the influence of doping on 
the AFM order. The size of the AFM moment changes continuously for lightly doped cases, 
but when the amount of doped carriers exceeds a certain value the AFM state is suddenly 
destroyed. We also show that Fermi surfaces remain and change significantly on doping even in 
the AFM state. This behaviour is explained by considering the nesting due to the multi-sheet 
Fermi-surface structure and multiorbital nature of the electronic bands characteristic to Fe 
pnictides. 

KEYWORDS: iron pnictides, magnetic order, orbital order, lattice distortion, Fermi surface, superconduc- 
tivity 



Since the discovery of superconductivity in 
LaFeAsOi-^F^ with a high transition temperature 
T c = 26 K, 1 ' extensive studies have been done on Fc 
pnictides. The main interests on these materials are 
not only on the high transition temperates such as 
T c = 55 K in SmFeAsOi_ x F x 2 ) and T c = 56 K in 
Gdi^^Th^FeAsO, 3 - 1 but also on the mechanism of the 
superconductivity. The electronic structure is quasi-two- 
dimensional 4-9 ) and superconductivity occurs around 
the magnetic phase boundaries 1 ' 10 ~ 13 ) as in bigh-T c 
cuprates. Such similarities suggest that magnetism is 
probably playing an important role in the emergence of 
superconductivity, and it is highly desirable to unveil 
the microscopic origin of magnetism characteristic to Fe 
pnictides. 

The magnetism in Fe pnictides is much different from 
that in cuprates. In the latter, ordering vector of the 
antiferromagnetism is (7r,7r), while it is (tt, 0) in Fc 
pnictides in the unfolded Brillouin zone (BZ) with one 
Fe ion per unit cell. 14-17 - 1 The undoped antiferromag- 
netic (AFM) states are metallic 10-12 ' 18 ~ 20 ) in Fe pnic- 
tides while insulating in cuprates. The AFM transition 
occurs a t 15 ~ 17 > 21 > 22 ) or near 14 ) the structural transition 
temperature in Fe pnictides. 

Such differences in magnetism may originate from 
the multiorbital electronic states and multi-sheet Fermi- 
surface (FS) structure in Fe pnictides. 4 ~ 9 ' Indeed, the 
AFM order with (tt, 0) due to nesting between hole and 
electron pockets [see Fig. 1(a)] has been suggested by 
using tight-binding models 23-25 ' and by band-structure 
calculations. 6 ' 7 ' Yildirim 26 ' has shown that the tetrago- 
nal lattice distortion occurs in the AFM state with (tt, 0), 
but not in the normal state. In addition, the lattice dis- 
tortion occurs neither in an AFM state with (tt, tt) nor 
in a ferromagnetic state. Ran et aZ. 27 ' have shown that 
a full band gap does not open and the system remains 
metallic even in the AFM state from a topological view 
point of the multiorbital system. 
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Fig. 1. (Color online) (a) Fermi surfaces of the two-orbital 
model 24 ) in the unfolded BZ at electron number n = 2 per site. 
The arrows indicate the nesting vector Q = (tt, 0). The hole sur- 
faces locate around (0,0) and (it, it), and the electron surfaces 
locate around (n, 0) and (0, tt). (b) Schematic views of the orbital 
states in the normal state (left) and the AFM state accompany- 
ing ferro-orbital order (right). The orbital stretched along x (y) 
axis represents the d zx (d yz ) orbital. The arrows represent the 
spin state, and the sizes of the orbitals indicate the occupancies 
of those orbitals. 



In this Letter, we show that such characteristic fea- 
tures on magnetism are explained even in the simplest 
model, i.e., by a two-orbital model, proposed for Fe 
pnictides 24 ' 25 -* by applying Hartree-Fock approximation. 
While the Hartree-Fock approximation has already been 
applied to the two-orbital model, 28 ' possibility of orbital 
order and effects of doping have not been investigated 
yet. They are main topics of the present study. The two- 
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orbital model cannot reproduce well the band structure 
obtained with the density functional theory, while a five- 
orbital model does well. 23 ) In addition, d xy orbital, which 
is not included in the two-orbital model, contributes 
some portions of the FSs. 29 ) Thus, to describe some prop- 
erties, all the five <i-orbitals may be necessary. However, 
the two-orbital model can reproduce at least the char- 
acteristic FS topology in Fe pnictides in the folded BZ, 
and is enough for the purpose of the present Letter. In 
this study, we take into consideration orbital order on 
an equal footing with AFM order, since they are closely 
related to each other in Fe pnictides. Indeed, the AFM 
order with the ordering vector Q — (n, 0) inevitably ac- 
companies ferro-orbital (FO) order [schematically shown 
in Fig. 1(b)] which results in a secondary orthorhom- 
bic distortion. In addition, we investigate doping effects 
on the antiferromagnetism and FSs reconstructed by the 
AFM order. As shown in Fig. 1(a), there are two kinds 
of nesting with the same nesting vector Q — (tt, 0), i.e., 
between FSs around (0,0) and (tt, 0), and between FSs 
around (0, it) and (tt,tt). The existence of the two kinds 
of nesting is important for stabilization of the AFM state 
against doping. We find that the structure of FSs changes 
significantly with doping even in the ordered state due 
to the multi-sheet FS nesting. 

In the two-orbital model, we consider a square lattice 
of Fe ions with d zx and d yz orbitals. 24 ' 25 ) The model 
Hamiltonian is given by 

h.r.r' ,a i,r 

+ U nj x njy + J cl xcr cl ya ,Ci xa 'Ci yi7 q\ 

+ J J]] C iTl C lTl CiT 'l CiT 'h 

where Ci ra is the annihilation operator of the electron 
at site i with orbital r and spin a (=| or j) and Ckra 
is the Fourier transform of Ci Ta . r = x and y represent 
d zx and d yz orbitals, respectively. rii Tcr = c\ TrT Ci Ta and 
rii T = n irtT- The coupling constants U, U', J, and J' 
denote the intraorbital Coulomb, interorbital Coulomb, 
exchange, and pair-hopping interactions, respectively. 
For the ti g orbitals, relations U = U' + J + J' and J = J' 
hold 30 ) and we use them. For the kinetic energy term, 
we use the hopping parameters proposed by Raghu e£ 
a/.: 24 ' tkxx = — 2£iCOsfc x — 2£2COsfc a — 4£3 cos k x cos k y , 
£kyy — —2t2Cosk x — 1t\ cos k y — 4^3 cos k x cos k y , and 
£kx V = £kyx = -4t4sinfc x sinfcy, where t\ = — 1,*2 = 
1.3t, £3 = ti = — 0.85£, and we have set the lattice con- 
stant unity. 

In this study, we consider weakly correlated cases, e.g., 
U/W ~ 0.29 for U/t = 3.5, where W = 12£ is the band- 
width. Thus, it is reasonable to apply a Hartree-Fock 
approximation. We assume that the expectation value of 
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Fig. 2. (Color online) (a) Order parameters as functions of U at 
n = 2 and J = 0.1(7. (b) Order parameters as functions of n at 
U/t = 3.5 and J = 0.1(7. 

the number ri{ ra is given by the following form: 
(»w) ={[n + m s (5 a i - S^i) + m (S TX - S Ty ) 
+ m so (S^ - S al )(S TX - 8 Ty )] 

(2) 

+ [n q + m sq (S a i - 6 al ) + m oq (5 TX - S Ty ) 

+ m soq {S a f - S al )(S TX - 8 Ty )]e iq ^}/^ 

where q = (tt, n) or (n, 0) = Q, Ti denotes the position 
of site i, and n is the number of electrons per site. The 
order parameters are m s , m , m so , n q , m sq , m oq , and 
m soq . We determine the lowest energy state among the 
solutions of the Hartree-Fock approximation. In Eq. (2), 
we consider the z component for the orbital state, i.e., 

m o = Tf^t.r.T'.ai^ira^T'^r'a) with a = Z, where T a is 

the Pauli matrix. We also considered order parameters 
with a = x and y, and we found that the z-component 
ordered state with q = Q = (tt, 0) always has lower 
energy than the other ordered states within parameters 
we investigate here. 

Figure 2(a) shows U dependence of the order param- 
eters m (for FO order), m s Q (for AFM order), and 
m so Q (for antiferro-spin-orbital order) at n = 2 and 
J = 0.1U. We find that the other order parameters are 
zero. As shown in the inset, m S Q jumps to a finite value 
at U/t ~ 2.97, and m and m SO Q also have jumps to fi- 
nite values at the same point while they are small and 
not visible on the scale of Fig. 2(a). Thus, the transition 
to the AFM state is of first order. 
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Figure 2(b) shows doping dependence of the order pa- 
rameters at U/t = 3.5 and J = 0.1U. We have cho- 
sen this value of U so as to the AFM state is destabi- 
lized by the doping of ~ 0.1 as in experimental observa- 
tions. 1 ' 10 ' 12 ' 13 ' The AFM moment changes continuously 
with doping at first, but suddenly disappears at n ~ 1.92 
and 2.1. We obtained small values of m s Q, which are 
much smaller than the saturation value n (for n < 2) 
or 4 — n (for n > 2), as in experimental observations 
( m sQ = 1 is corresponding to 1 /Ub of an AFM moment): 
0.25 fin, 0.35 (i B , or 0.36 n B in LaFeAsO; 13 ' 14 ' 31 ' 32 ' 

0. 94 fi B or 1.01 fi B in SrFe 2 As 2 ; 15 ' 17) 0.4 fi B or 0.87 fi B 
in BaFe2As2. 16, 18 - > The estimated values of the ordered 
moments depend on the experimental probes even for the 
same material, and the reason of this discrepancy is not 
clear at present. In the electron doped case, the AFM 
moment increases with doping, and the transition tem- 
perature is expected to become higher than the undoped 
case. It is in contradiction to experimental observations, 
and to resolve this discrepancy, we have to extend the 
model, e.g., by using the five-orbital basis. 

Note that in the AFM state with Q = (tt,0), x and y 
directions are not equivalent, and the occupancies of d zx 
and d yz become different. Thus, the AFM state in the 
multiorbital system inevitably accompanies FO order, 

1. e., finite m . Through an electron-lattice interaction, 
the FO order results in a lattice distortion from a tetrag- 
onal to orthorhombic structure. This is consistent with 
experimental observations that the AFM phase is always 
orthorhombic. The obtained small values of the order pa- 
rameter m for the FO order may be responsible for the 
weakness of anomaly in lattice distortion, e.g., small vol- 
ume change in SrFe2As2 at the transition. 21 ' Note that 
another scenario is proposed for the lattice distortion, 
in which the lattice distortion relaxes magnetic frustra- 
tion and is necessary for occurrence of the AFM order. 26 ' 
On the contrary, in our theory, the lattice distortion is 
a secondary effect due to the FO order accompanied by 
the AFM order. Since the AFM state in Fe pnictides 
is metallic, we believe that our picture is more suitable 
for Fe pnictides. In the coexistent state of antiferromag- 
netism and FO order, m so Q also becomes finite as shown 
in Fig 2. Note that we obtain similar results for J = 
at least in a small- 1/ region and the choice of the value 
of J does not change the present results qualitatively. In 
this model, at zero temperature, the AFM order and the 
FO order disappear at the same doping and it is consis- 
tent with the experimental observations. However, some 
extensions, e.g., inclusion of electron-lattice interaction, 
may be necessary to obtain different transition temper- 
atures for the AFM order and for the FO order as in 
LaFeAsO. 

To obtain further insights into the ordered states, we 
show FSs in Fig. 3 in the normal and ordered states at 
half-filling (n — 2) and at around phase boundaries. It 
is evident from the figure that the x and y directions 
are not equivalent in the ordered states. In the hole 
doped case, n = 1.93, nesting between FSs centered at 
k = (0, 0) and (n, 0) in the unfolded BZ is strong. Then, 
these FSs are reconstructed into small pockets around 
(~ ±7r/4,0), while the other FSs centered at (0,7r) and 
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Fig. 3. (Color online) FSs in the folded BZ with respect to Q = 
(tt, 0) for n = 1.93, 2, and 2.07. The upper panels show the FSs 
in the normal state. The lower panels show those in the AFM 
state with FO order at U/t = 3.5 and J = 0.1(7. 



(•7T, tt) are almost unchanged. On the other hand, in the 
electron doped case, n — 2.07, nesting is strong between 
FSs centered at k = (0, tt) and (tt, tt). Around zero dop- 
ing, both types of nesting can contribute to stabilize the 
AFM state. As a result, the doping effect on the AFM 
moment is not so significant in the lightly doped cases, 
while the structure of the FSs changes very much. Thus, 
the multi-sheet FS nesting is important for the stabiliza- 
tion of the AFM state in this system. When we dope 
carriers further, the nesting becomes weak and the AFM 
state is destabilized suddenly. 

Note that such a mechanism to stabilize the AFM state 
against doping is applicable as long as the sizes of the 
hole surfaces around (0, 0) and (tt, tt) are different and 
the undoped system is a compensated or nearly com- 
pensated metal. The smaller hole surface mainly con- 
tributes for the realization of antiferromagnetism for the 
hole doped case, and the larger hole surface mainly con- 
tributes for the electron doped case. Thus, this mech- 
anism works irrespective of precise choice of the model 
parameters. We also note that in a more realistic model, 
i.e., five-orbital model, 23 ' both the two hole surfaces lo- 
cate around (0, 0) even in the unfolded BZ. However, the 
nesting vector Q = (tt, 0) is the same as in the two-orbital 
model, and the present mechanism to stabilize the AFM 
state is applicable provided the two hole surfaces have 
different sizes. 

At n = 2 only small pockets of FSs remain in the or- 
dered states. The area of one pocket at n — 2 is 0.86% of 
the folded BZ for Q = (tt, 0) [and of the normal state BZ 
folded due to the actual lattice structures of Fe pnictides 
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Fig. 4. Band structure in the folded BZ for Q = (tt, 0) and den- 
sity of states at n = 2. (a) Band structure and (b) density of 
states in the normal state, (c) Band structure and (d) density 
of states in the AFM state with FO order at U/t = 3.5 and 
J = 0.1(7. The Fermi energy is set to be zero in these figures. 



(two Fe ions per unit cell)]. There are two electron pock- 
ets and two hole pockets, but the areas of them are the 
same, since the model is a compensated metal at n — 2 
and the two electron (hole) pockets occupy the same 
amount of area due to symmetry. Experimentally ob- 
served volumes of FSs in the AFM state are small: 0.26%- 
1.38% in SrFe 2 As 2 33 ) and 0.3%-1.7% in BaFe 2 As 2 34 ) of 
the normal state folded BZ. These values are comparable 
with our theoretical ones. In the normal state at n = 2, 
the hole pocket around (0, 0), hole pocket around (n, tt), 
and electron pocket around (tt, 0) occupy 7.13%, 13.24%, 
and 10.18% of the folded BZ, respectively. These values 
are also comparable with experimental ones, 2.8%-9% of 
BZ in LaFePO in the normal state. 35 ) 

Figure 4 shows the band structure and density of states 
in the normal and ordered states. In the ordered state, 
band gaps open at some points at the Fermi level, while 
not at (~ 7r/4, 0) and (~ 7r/4, it). At k y = and tt, the 
off-diagonal element €k xy in the kinetic energy term is 
zero, and d zx and d yz orbitals do not mix. The mean field 
in the ordered state mixes electrons with k and k + Q 
in the same orbitals, and d zx - and d yz -orbital states are 
not mixed even in the ordered state at k y — and tt. 



The two bands crossing at around the Fermi level are 
different orbitals at both (~ 7r/4, 0) and (~ 7r/4,7r), and 
a gap cannot open there. As a result, the FSs do not 
disappear on the lines k y — and tt even in the AFM 
state as shown in Fig. 3. Note that in Fig. 3 the Fermi 
pockets at n = 2.07 on k y = tt are very small in the 
ordered state but have finite volumes. Thus, the system 
remains metallic in the AFM state as in experimental 
observations. The density of states in the ordered state 
has a gap-like structure around the Fermi level, but re- 
mains finite at the Fermi level. For larger U cases, the 
band structure changes very much, and the system can 
become insulating. 28 -* By increasing U, the band crossing 
points on (0, 0)-(n/2, 0) and (0, tt)-(tt/2, tt) lines move to 
(7r/2,0) and (7r/2,7r), respectively, and finally the band 



crossing disappears. Then, the system can become insu- 
lating, e.g., at U/t > 6.7 for n = 2 and J = 0.1U. Thus, 
U should not be very large in Fe pnictides. 

In conclusion, we have shown that characteristic fea- 
tures of the AFM state in Fe pnictides can be natu- 
rally understood within the two-orbital model. The sta- 
bility of AFM phase is due to the multi-sheet FS nest- 
ing. The tetragonal to orthorhombic lattice distortion 
is a secondary effect due to the FO order but not a 
driving mechanism of antifcrromagnctism. FSs remain 
in the ordered state due to the multiorbital character of 
the crossing bands. Our theory indicates that the FSs 
change significantly upon doping. In the doped AFM 
states around phase boundaries, some Fermi pockets be- 
come very small, while the other Fermi pockets have large 
volumes as in the normal state. Experimental observa- 
tions of these FSs are highly desired, since we can know 
what kind of nesting is strong around the phase bound- 
aries from the reconstructed FSs. Such a knowledge is 
important to unveil fluctuations which mediate the su- 
perconducting pairing. 
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